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Abstract 

A search for resonances decaying to top quark-antiquark pairs is performed using a 
dilepton+jets data sample recorded by the CMS experiment at the LHC in pp colli- 
sions at y/s = 7TeV corresponding to an integrated luminosity of 5.0fb _1 . No sig- 
nificant deviations from the standard model background are observed. Upper limits 
are presented for the production cross section times branching fraction of top quark- 
antiquark resonances for masses from 750 to 3000 GeV. In particular, the existence of a 
leptophobic topcolor particle Z' is excluded at the 95% confidence level for resonance 
masses M z < 1.3 TeV for T z , = 0.012M Z /, and M z > < 1.9 TeV for T z , = 0.10M Z /. 
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1 Introduction 

Electroweak symmetry breaking is a cornerstone for the understanding of particle physics. 
However, despite the spectacular phenomenological success of the standard model (SM), and 
the recent observation of a new boson at the Large Hadron Collider (LHC) [1, 2j, the precise 
mechanism of electroweak symmetry breaking remains unknown. Various new models have 
been proposed to explain this mechanism. One such class of models, topcolor-assisted techni- 
color (TC2) |3]-(5|, provides a dynamical explanation for electroweak symmetry breaking and 
flavor symmetry breaking, giving masses to the weak gauge bosons and fermions. Under one 
of the scenarios of TC2, a heavy boson Z' is predicted with preferential couplings to the third 
quark generation and with no significant couplings to the leptons ("leptophobic"). 

Direct searches for massive resonances that decay preferentially to top quark-antiquark (tf) 
pairs are currently feasible only at hadron colliders. Experiments seek to observe an excess be- 
yond that predicted by the SM, typically in the distribution of the invariant mass of the tt decay 
products. Searches in pp collisions at the Tevatron and the early searches in pp collisions at the 
LHC by the ATLAS experiment have excluded a narrow-width, leptophobic Z' with a mass 
lower than 900 GeV HrEJ. The searches by the Compact Muon Solenoid (CMS) experiment at 
the LHC have excluded a narrow-width, leptophobic Z' in the mass range 1.3-1.5 TeV and in a 
narrow window around 1 TeV |9lH0l|. This paper describes a search for a Z' — > tt resonance in 
pp collisions in the 2.1 + 2v + jets final state, where I is an electron (e) or a muon This is the 
first search for topcolor leptophobic Z' in final states involving two leptons. The data sample 
corresponds to a total integrated luminosity of 5.0 fb _1 [11 J at y/s = 7 TeV collected by the CMS 
detector in 2011. 

2 The CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid, 13 m in length and 
6 m in diameter, which provides an axial magnetic field of 3.8 T. The bore of the solenoid is 
surrounded by various particle detection systems. Charged particle trajectories are measured 
by a silicon pixel and strip tracker, covering < cp < In in azimuth and \t]\ < 2.5, where the 
pseudorapidity rj is defined as rj = — ln[tan(0/2)], and 9 is the polar angle of the trajectory 
of the particle with respect to the counterclockwise-beam direction. A crystal electromagnetic 
calorimeter and a brass/ scintillator hadronic calorimeter surround the tracking volume. The 
calorimeter provides high-resolution energy measurement of electrons. Muons are measured 
in gas-ionization detectors embedded in the steel flux return yoke outside the solenoid. The 
detector is nearly hermetic, which facilitates the measurement of energy balance in the plane 
transverse to the beam direction. A two-tier trigger system selects the most interesting pp 
collision events for use in physics analysis. A more detailed description of the CMS detector 
can be found in Ref. |T2|. 

3 Event reconstruction 

In the Z' — > tt search, a tf decay topology is used where each top quark decays to a W boson 
and a b quark, and subsequently each W boson decays into a lepton and a neutrino. The 
signature for such an event is two oppositely charged, isolated leptons with high transverse 
momenta (px), large momentum imbalance due to two undetected neutrinos, and at least two 
jets. Events are required to pass a trigger requiring at least two high-px isolated leptons, and are 
separated into three channels based on lepton flavor: ee, and efi. The principal sources of 
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4 Signal and background modeling 



background are SM tt, Z/7* -4 (Drell-Yan, DY), single-top-quark, and diboson (WW, WZ, 
and ZZ) production. Other minor contributions are from W — > Iv and multijet production. 
Electrons, muons, jets, and the momentum imbalance are reconstructed using a particle-flow 
algorithm [ 13 1 . The negative of the vector sum of the momenta of all reconstructed particles in 
the plane transverse to the beams is the missing transverse momentum [Tl4| , whose magnitude 
is called missing transverse energy (^t). The identification criteria of each object and additional 
selections are chosen to reduce all backgrounds other than SM tt production. 

Electron candidates are reconstructed from clusters of energy deposits in the electromagnetic 
calorimeter, which are then matched to hits in the silicon tracker. Electron identification is 
based on shower-shape and track-cluster matching variables ffl5| . Electrons are required to 
have px > 20GeV and \tj\ < 2.5, and are excluded if they are in the transition region between 
the barrel and endcap calorimeters, 1.4442 < < 1.5560, because their reconstruction in this 
region is degraded due to additional material there. The electron track must be within 0.04 cm 
of the primary vertex in the plane transverse to the beam. Additionally, electrons coming from 
photon conversions in the detector material are rejected if there are missing hits in the inner 
tracker layers or if there is another close track with opposite charge and with a similar polar 
angle. 

Muons are reconstructed using the information from the muon detectors and the silicon tracker [TBI . 
The reconstructed muon track must be within 0.02 cm of the primary vertex in the plane trans- 
verse to the beam. Muons are required to have pj > 20 GeV and \rj\ < 2.4. 

To remove leptons arising from decays of hadrons immersed in jets, the electrons and muons 
are required to be isolated. The isolation requirement is based on the ratio of the total transverse 
energy observed from all particles in a cone of size R = \J (A(p) 2 + (A//) 2 < 0.3 centered on the 
direction of the lepton to the transverse momentum of the lepton. This quantity must be less 
than 0.17 for electrons, and less than 0.20 for muons. 

In order to reduce the background from low-mass dilepton resonances, events are rejected if the 
dilepton invariant mass Mm < 12 GeV. To suppress the contribution from Z-boson production, 
a veto on events with 76 < Mm < 106 GeV is applied in the ee and fifi channels. 

Events are required to contain at least two jets, reconstructed using an anti-Zcj clustering algo- 
rithm with a distance parameter of 0.5 |T7). Corrections are applied to account for the depen- 
dence of the detector response to jets as a function of 77 and pj and the effect of pileup (multiple 
pp collisions) [18J. The corrections are propagated to recalculate the missing transverse energy. 
The jets are required to have pj > 30 GeV and \rj | < 2.5. Additionally, at least one of the jets is 
required to be tagged as a b jet based on the identification of a secondary vertex |Tl9| . 

Finally, in order to further reduce the DY and multijet backgrounds in the ee channel, and the 
DY contribution in the channel, a requirement of #j > 30 GeV is applied in these channels. 
The DY and multijet backgrounds are negligible in the eji channel. 

4 Signal and background modeling 

The signal efficiency and background rejection of the selection outlined above are determined 
from simulation studies augmented where necessary by corrections based on control samples 
in data. The resonance signal Z' — > tt is modeled using the MadGraph 5.1.1 pOfl Monte 
Carlo (MC) event generator, with the top-quark mass set to 172.5 GeV and CTEQ6L |2"T| parton 
distribution functions. Samples are generated with resonance masses between 750-3000 GeV, 
and for two resonance-width scenarios: narrow (Tz' = 0.012Mz') and wide (Tz> = O.lOMz')- To 
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calculate the expected number of signal events, we use cross sections for a leptophobic topcolor 
71 p2|. A scale factor of 1.3 is used to account for next- to-leading-order (NLO) corrections p3|. 

The background events from SM tt, DY, and W — >• Iv are generated using MadGraph 5.1.1. 
Diboson events are generated using pythia 6.424 fl24| . and single-top-quark events are gen- 
erated using POWHEG 1.0 rl380 [25-27]. The MadGraph and POWHEG events are processed 
through PYTHIA in order to add the initial- and final-state radiation and showering, together 
with the production of the underlying event [28]. To estimate the expected number of back- 
ground events, the background samples are normalized to the theoretical cross sections shown 
in Table [l] All MC events are processed through a simulation of the CMS detector based 

Table 1: Theoretical cross sections, including higher-order corrections, for the SM back- 
grounds Il29ti33l. 



Background 


Cross section [pb] 


tt 


160 


DY (M(f > lOGeV) 


15000 


WW 2£2v 


4.5 


WZ -4 3£v 


0.61 


ZZ (inclusive) 


7.4 


W 


31000 


Single top 


85 



on GEANT4 [34], and are overlaid with events from minimum-bias interaction to account for 
pileup effects at high instantaneous luminosity. An additional set of corrections is applied to 
account for data-MC differences. These include reweighting MC events to match the overlaid 
pileup distribution to that inferred from the data, and applying a scale factor of 0.95 per tagged 
b jet [19] to account for the observed difference in b-tagging efficiency between data and MC. 

The simulation of the DY background does not adequately reproduce the production rate, es- 
pecially in the presence of missing transverse energy. The overall normalization for the DY 
process is therefore obtained using data in the DY-enriched region of 76 < < 106 GeV in 
the ee and channels. This region is excluded from the dataset used for the main analysis 
and thus provides an independent DY control sample. In the eji channel, the DY-enriched re- 
gion is between 40 < Ma < 70 GeV since the main contribution to DY in this channel is from 
Z/ 7* — > tt events. These events have a peak at slightly lower dilepton masses on account 
of the invisible decay products of the t lepton. As the Z-mass veto is not applied in the eji 
channel, the requirement of at least two jets is modified to exactly one jet in order to ensure the 
exclusion of this calibration sample from the main signal sample. The normalization factors 
obtained are 1.34, 1.24, and 1.20 in the ee, and efi channels, respectively. 

The multijet background is estimated directly from data. This background is from misidentified 
leptons or genuine leptons from semi-leptonic decays of b/b or c/c quarks, which pass the 
isolation requirement. It is determined from data by inverting the isolation criteria for both 
leptons, and then extrapolating that yield to the signal region. 



5 Event yields 

The number of events for the expected SM backgrounds and the observed data after all selec- 
tions and applying corrections to account for data-MC differences are listed in Table [2] The 
uncertainties shown for the backgrounds are from the systematic effects discussed in Section|7| 
There is good agreement between data and SM backgrounds in all three channels. For com- 
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parison, the number of expected 71 signal events for Mz> = 750-3000 GeV are listed in Table [3] 
Table 2: Event yields for SM backgrounds and data in the ee, and ep channels. 



Sample 


ee 


W 


ep 


tt 


2210 ± 460 


2550 ± 550 


7300 ± 1500 


DY 


410 ± 130 


420 ± 130 


179 ± 57 


Diboson 


11.5 ± 1.5 


15.4 ± 2.0 


32.3 ± 4.2 


W ->■ Iv 


17 +24 


q+3.4 


zo -26 


Single top 


106 ± 15 


121 ± 18 


343 ± 50 


Multijets 


41.8 ± 6.9 


50 ± 10 


103 ± 14 


Total background 


2790 ± 510 


3150 ± 600 


8000 ± 1500 


Data 


2690 


3098 


7704 



Table 3: Event yields for a leptophobic 71 in ee, pp, and ep channels. 



Sample 


ee 


W 


ep 


T z >/M z > 


= 0.012 




71 (750 GeV) 


67 


79 


200 


71 (1000 GeV) 


26 


28 


68 


71 (1250 GeV) 


8.2 


9.8 


22 


71 (1500 GeV) 


2.9 


3.1 


7.0 


71 (2000 GeV) 


0.3 


0.4 


0.8 


r Z '/M z 


/ = 


0.10 




71 (1000 GeV) 


180 


200 


480 


71 (1500 GeV) 


23 


26 


57 


71 (2000 GeV) 


2.9 


2.9 


7.0 


71 (3000 GeV) 


0.1 


0.1 


0.2 



Distributions of the transverse momentum of the highest-px electron in the ee channel, the 
pseudorapidity of the highest-px jet in the pp channel, and A(p between the missing transverse 
momentum and highest-px lepton in the ep channel are shown in Fig. [I] There is good agree- 
ment between data and the sum of all SM backgrounds. Similarly, a good agreement is seen in 
all four-vector distributions (px, rj, and <p) of all final-state objects in the three channels. 

The tt invariant mass is constructed using the four-vectors of the two leading leptons, the two 
leading jets, and the missing transverse energy. The longitudinal momenta p 2 of the two neu- 
trinos in the final state cannot be measured experimentally and are set to zero. The tf invariant 
mass distributions for the data, the sum of all backgrounds, and the 71 signals for the narrow- 
width scenario (T Z / = 0.012Mz/) are shown in Fig. [2]for the ee, pp, and ep channels. The data 
is described well by the SM backgrounds, and there is no statistically significant evidence for 
the presence of a 71 signal. 

6 Bayesian neural network analysis 

A multivariate analysis, based on Bayesian neural networks (BNN) [35 1, has been carried out to 
provide a more powerful discriminant between backgrounds and the 71 signal than that based 
on invariant mass alone. The discriminant lies in the interval [0,1]. It is constructed such that 
the signal events tend to have a value closer to 1 than the background events that peak instead 
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Figure 1: Distributions of the transverse momentum of the highest-px electron in the ee chan- 
nel, the pseudorapidity of the highest-px jet in the pp channel, and Acp between the missing 
transverse momentum and highest-px lepton in the ep channel. The hatched region indicates 
systematic uncertainties on the sum of SM backgrounds. 




M(tt) [GeV] M(tf) [GeV] M(tt) [GeV] 

Figure 2: Distributions of the tt invariant mass for the ee, pp, and ep channels. The p z values 
for both neutrinos are set to zero. The hatched region indicates systematic uncertainties on the 
sum of SM backgrounds. The Z' signal corresponds to a resonance width of T z > = 0.012Mzv 
and has been scaled up so as to be visible. 

closer to 0. To build this discriminant the background is defined as the sum of all SM processes, 
and the signal is set to hAjj = 750 GeV for narrow width, and Mz> = 1000 GeV for wide width. 
The separation power increases with hypothetical Z' mass. Thus, using the Z' simulated sample 
with the lowest Z' mass for the signal when training the BNN ensures good discrimination 
between signal and background even for the higher Z' masses. A separate discriminant is 
constructed in each of the three channels. As inputs to the training, the following 17 variables 
are used in each of the two scenarios, and in each channel: 

• px and rj of the highest-px lepton, 

• px, f], and Acp of the second highest-px lepton, 

• px, f], and Acp of the highest-px jet, 

• px, f], and A(p of the second highest-px jet, 

• Ej, and Acp of the missing transverse momentum, 

• px, J], and Acp of the highest-px b-tagged jet, and number of b-tagged jets, 

where Acp is the difference in azimuth between the object and the highest-px lepton. This set of 
input variables constitutes the full array of four vectors of final-state objects that are measured 
in the analysis, along with additional information about the b-tagged jets. 



6 



7 Systematic uncertainties 



The resulting BNN outputs for the observed data, the SM background, and the Z' signals for 
Tz' = 0.012Mz' and Tz> = O.lOMz' are shown in Figs. [3] and |3J respectively, for the ee, fifi, 
and efi channels. There is good agreement between data and the SM background in all three 
channels with no evidence of a resonance signal. Upper limits are set on cr Z 'B(Z' — > tt). 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

BNN output BNN output BNN output 



Figure 3: Distributions of the BNN output discriminant for the ee, fiji, and efi channels. The 
hatched region indicates systematic uncertainties on the sum of SM backgrounds. The Z' signal 
corresponds to a resonance width of T Z / = 0.012Mz', and has been scaled up so as to be visible. 
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BNN output BNN output BNN output 



Figure 4: Distributions of the BNN output discriminant for the ee, fifi, and efi channels. The 
hatched region indicates systematic uncertainties on the sum of SM backgrounds. The Z' signal 
corresponds to a resonance width of T Z ' = O.lOMz/, and has been scaled up so as to be visible. 



7 Systematic uncertainties 

The signal and background models are affected by a number of systematic uncertainties, which 
are propagated into the limit calculation. The uncertainties are divided into two categories: 
those that affect only the overall normalization of a process ("rate") and those that affect also 
the distribution of the BNN discriminant ("shape"). The rate effects include the uncertainty on 
predicted cross section and normalization for each SM background based on data, as discussed 
in Section |4j and the uncertainties from integrated luminosity, lepton identification and isola- 
tion, b-tagging scale factor, jet energy scale, and pileup reweighting for both SM background 
and Z' signals. Rate uncertainties are also included for the SM tt and W — > iv events from 
variations in the renormalization and factorization scales (A), and the matching scale for jet 
production threshold between jets from matrix-element generation in MADGRAPH and parton 
showering in PYTHIA [36 J. The A and matching scales are each varied up and down by a factor 
of two with respect to their nominal values (Q 2 and 20GeV, respectively) in order to estimate 
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the uncertainty. The shape effects include the change in shape of the BNN distributions from 
the uncertainties from jet energy scale, pileup reweighting, and A and matching scales. The 
uncertainty due to parton distribution functions is negligible and therefore not included. All 
uncertainties are summarized in Table |4j where each row represents an independent entity. The 
total uncertainty on the sum of all SM backgrounds is 18%. 

Table 4: Rate (top part) and shape (bottom part) uncertainties, from different components, 
affecting the total expected number of signal and background events. The uncertainties on the 
A and matching scales in SM tt, the jet energy scale, and the pileup reweighting are indicated 
by their range across bins of the BNN distribution. 



Component Uncertainty [%] 

tt cross section 15 
DY normalization 30 
Diboson cross section 3.8 
W — > Iv cross section 5.0 
Single top cross section 7.7 
QCD normalization 13 (ee), 18.2 (jiji), 9.7 (e^) 

Integrated luminosity 2.2 
Lepton identification 2.0 
b-tagging scale factor 10 
A scale (W -> iv) 100 
Matching scale (W Iv) 100 

A scale (tt) 1.9-2.9 

Matching scale (tt) 3.4-5.5 

Jet energy scale 0.3-4.7 

Pileup reweighting 0.3-2.4 



8 Results 

With no excess observed, upper limits on az'B(Z' — > tt) at the 95% confidence level (CL) for 
different values of M Z / are set using the CL S criteria H37ll38| . All systematic effects are included 
with correlations across the different samples and channels. The sensitivity of the results is 
estimated using the invariant mass distributions shown in Fig.[2]and the BNN output distribu- 
tions shown in Fig.[3j and comparing the expected limit on djjB^Zl — > tt) for the two methods 
at Mz' = 750 GeV and r Z / = 0.012M Z /. An expected limit is obtained from an ensemble of sim- 
ulated pseudo-data sets, where each set is constructed from the background-only hypothesis. 
Using the BNN distribution improves the expected limit by 29% compared to using invariant 
mass distribution. For this reason, the more sensitive BNN technique is used for the subsequent 
measurements. The resulting expected limits and the observed limits using data are shown in 
Fig.|5]for both narrow and wide resonances. 

The theoretical predictions for a leptophobic Z' p2ll23l are used to exclude heavy Z' resonances 
of masses M z > < 1.3 TeV for a width T z , = 0.012M Z <, and M Z / < 1.9 TeV for a width T z , = 
0.10M Z /. In the current analysis, the expected lower limits on M z > are 1.1 TeV and 1.7TeV for 
r z , = 0.012M Z / and T z , = 0.10M ZV respectively. 
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Figure 5: Upper limits on a Z iB{Z! — > tt) versus M Z i for narrow and wide resonances. Also 
shown is the theoretical prediction for the Z'. 



9 Summary 

A data sample, corresponding to an integrated luminosity of 5.0 fb -1 collected in pp collisions 
at a/s = 7 TeV, has been analyzed in a search for heavy resonances decaying to top quark- 
antiquark pairs with subsequent leptonic decay of both top quark and antiquark. No excess 
beyond the standard model prediction is observed. Upper limits at the 95% CL are derived on 
the product of the production cross section and branching fraction for these decays, for various 
masses of narrow and wide resonances. The existence of a leptophobic 2! topcolor particle is 
excluded for M z , < 1.3 TeV with T z , = 0.012M ZV and for M z , < 1.9 TeV with T z , = O.WM z >. 
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